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We investigate the sedimentation equilibrium of a charge stabilized 
colloidal suspension in the regime of low ionic strength. We analyze the 
asymptotic behaviour of the density profiles on the basis of a simple 
Poisson-Boltzmann theory and show that the effective mass we can 
deduce from the barometric law corresponds to the actual mass of the 
colloidal particles, contrary to previous studies. 



I. INTRODUCTION 

Under the action of gravity a colloidal suspension sediments to form a stratified 
fluid. The equilibrium density profile of the colloidal particles results from the 
balance between the gravitational force and thermodynamic forces as derived from 
the free energy of the system. The density profiles usually exhibits a dense layer 
of colloidal particles at the bottom of the container above which a light cloud of 
colloidal particles floats. In this last regime, the density of particles is small enough 
to treat the fluid as an ideal gas. Under the reasonable assumption that density 
gradients can be neglected, the equilibrium colloidal density obey the well known 
barometric law: 

pcolO) = p° col exp(-z/l g ) (1.1) 

Here, p co i(z) denotes the density profile of the colloidal particles, z is the altitude 
and l g — (PMg)^ 1 is the gravitational length where (3 — (fceT) -1 is the inverse 
temperature, M is the buoyant mass of a colloidal particle and g the intensity of 
the gravitational field. This exponential law is of practical interest since it gives 
a prescription for the measurement of the buoyant mass M of the particles. How- 
ever a recent experimental study of the sedimentation profiles of strongly de-ionized 
charged colloidal suspensions jjj lead the authors to challenge the validity of this 
barometric law. An exponential behaviour was indeed observed in the asymptotic 
regime, but the measured gravitational length I* could differ significantly from the 
expected one (a factor of two). I* was found to systematically overestimate the 
actual value l g , with the result that the buoyant mass measured within these ex- 
periments is systematically reduced compared to the known buoyant mass of the 
particles. 

Some theoretical efforts have been made to study this problem. First Biben and 
Hansen p) solved numerically the problem in a mean field approach, but unfor- 
tunately due to numerical difficulties the samples height considered where of the 
order of the micron while in the experiment the samples height are of the order of 
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the centimeter. As a consequence, the dilute region at high altitude could not be 
studied in this approach. Nevertheless the numerical results show a positive charge 
density at the bottom of the container and a negative charge at the top while the 
bulk of the container is neutral. This result show that a non-zero electric field exists 
in the bulk of the container and acts against gravity for the colloids. 

More recently one of the authors studied a two-dimensional solvable model for 
this problem ||. This model is not very realistic (the valency of the colloids was 
Z = 1 and there was no added salt) but has the nice feature of being exactly solvable 
analytically. It confirmed the condenser effect noticed for small height containers 
in Ref. ||]. For large height containers it showed a new interesting phenomenon: 
while there is still a positive charge density at the bottom of the container, the 
negative charge density is not any more at the top of the container floating but 
at some altitude. Interestingly, the analytical expression for the density profiles in 
the asymptotic regime predicts a decay in exp(— z/l g )/z for the colloidal density. 
Besides the 1/z factor that cannot be explained by a mean field approach, no mass 
reduction is predicted by this model. However one should be cautious when com- 
paring two-dimensional systems to the three dimensional case because the density 
in not relevant in two-dimensional Coulomb systems: no matter how small the den- 
sity is the system is always coupled, the ideal gas regime is never attained. For this 
reason a decay of the density similar to the one of an ideal gas is in itself surprising 
in two dimensions. 

Lately new results based on an approximate version of the model introduced in 
reference Q lead the authors of these studies to conclude that the mean-field 
approach was indeed able to predict a mass reduction in the asymptotic regime. 
Here we present some new results about this problem treated under the Poisson- 
Boltzmann approximation, and show that it is indeed not the case. 

II. THE MODEL AND THE POISSON-BOLTZMANN APPROXIMATION 

Let us consider some colloidal particles (for example some latex spheres) in a 
solution with some amount of added salt. In a polar solvent like water the colloids 
release some counterions and therefore acquire a surface electric charge Ze {Z is a 
entire number usually positive and — e is the charge of the electron). We consider 
that the colloidal sample is monodisperse, all colloids have the same valency Z, and 
that the counterions and the salt cations are both monovalent and therefore we shall 
not make any distinction between cations coming from the colloids and salt cations. 
We then consider a three-component system composed of colloidal particles with 
electric charge Ze and mass M, counterions with charge — e and coions with charge 
+e. We shall neglect the masses of the counterions and coions when compared 
with the mass of the colloids. The solvent shall be considered in a primitive model 
representation as a continuous medium of relative dielectric permittivity e (for water 
at room temperature e rs 80). The system is in a container of height h, the bottom 
of the container is at z = altitude. We consider that the system is invariant 
in the horizontal directions. The density profiles of each species are denoted by 
p co \(z), p+(z) and p-(z) (z is the vertical coordinate) for the colloids, the cations 
and the anions respectively at equilibrium. Let us define the electric charge density 
(in units of e) p = Zp co \ — p_ + p + and the electric potential <f>, solution of the 
Poisson equation 

d 2 $ , , Am , . . . 

-^{z) = -—ep{z) (2.1) 

It is instructive to recall that the Poisson-Boltzmann equation can be derived 
from the minimization of the free energy density functional 
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f / Mgzp col {z)dz + - J ep(z)$(z) (2.2) 



where is the de Broglie wavelength of species i. Minimization of the grand 
potential with respect to the densities: SJ- '/ 'Spi(z) — pi = 0, where is the chemical 
potential of species i, yields 

p co i(z) = p° col exp(-/3Ze$(z) - (3Mgz) (2.3a) 
p+(z) = p° + cxp(-/3e$(z)) (2.3b) 
p- (z) = p°_ exp(/?e$(z)) (2.3c) 

We shall work in the canonical ensemble, the prefactors p° which depend on the 
chemical potentials pn are determined by the normalizing conditions 

h 

p i {z)dz = N i (2.4) 







where Ni is the total number of particles per unit area of species i. The system is 
globally neutral so we have ZN co \ — iV_ + N + = 0. 

Let us introduce the following notations: l g — (fiMgj^ 1 is the gravitational length 
of the particles, I = (3e 2 /e is the Bjerrum length, (j> = {3e$ is the dimensionless 
electric potential and Ki = (^TrlNi/h) 1 ^ 2 . k^ 1 are the Debye lengths associated to 
the counterions and the coions and (Zk co i) _1 is the Debye length associated to the 
colloidal particles. For a quantity q(z) depending on the altitude, let us define its 



mean value (q) — J Q q(z) dz/h. With these notations equations (2.1) and (2.2) yield 
the modified Poisson-Boltzmann equation 

d^ r , _ _ 3 e-f*W-^_ a e*W a e-+W 



From Eq. (2.5) it is clea r that the problem has the following scale invariance: if 
<j>(z) is a solution of (|2.5| ) then (f>(az) is a solution of the problem with the rescaled 
lengths al g and an^ . 

The advantage of the density functional formulation of the problem is that it 
allows for systematic corrections to the Poisson-Boltzmann approximation. For 
instance, one may be interested in the effect of the finite size of the macroions. Let 
a be the diameter of the colloids, 7y co i = 7ra 3 /9 co i/6 the volume fraction of colloids 
and p± = p±/ (1 — ?7 C oi) the effective densities of the microions. Then, the finite size 
of the colloids can be accounted in a local density approximation (LDA) by adding 



to the free energy density functional (2.2) the free energy excess term 

/cxcc(p C oi(z)) dz (2.6) 

where / eX ec(pcoi) is the excess free energy of a hard sphere fluid derived by the 
Carnahan-Starling equation of state 

/cxcc(Pcol) = fcBrPco! 4 ^ 01 " 377 ,! 1 (2.7) 



and by replacing in (2.2) the entropy term of the microions by 
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k B T / p±{z){hx(\ 6 ± p* ± {z)) -l)dz. 



Mini miza tion of this new free energy functional gives the modified version of 
Eqs. (U) 



PcoiO) = p° col exp{-f3Ze$(z) - f3Mgz) 

"3?7 col (z) 3 - 9t] co1 (z) 2 + 8?7 col (z) ttct 3 



x exp 



+ — (p* + (z) + p*_(z)) 



(2.9a) 



(1 - Vcoi(z)) 3 6 
p+(z) = p° (1 - Vcol (z j) cxp(-/?e$(z)) (2.9b) 
p_(z) =p°_(l-7? col (z))exp(/3e$(z)) (2.9c) 

There are several length scales in this problem: the gravitational length of the 
colloids l g , the Debye or screening length, the height h of the container and eventu- 
ally the hard core diameter of the particules. In a realistic case h is of the order of 
the centimeter, l g of the order of 0.1 mm, the screening length of the order of 10 nm. 
We are faced to a practical numerical problem, when we will transpose the problem 
to a lattice, the lattice spacing should be smaller than all the physical lengths, but 
since h is much larger than the others lengths, the number of sites in the lattice 
should be very high (of order 10 6 ). A possible approach to deal with this problem 
is to study small containers as in Ref. @. In this paper we want to study the case 
of high containers so we will consider very deionized systems in which the screening 
length is of the order of 0.1 mm, much larger than in usual physical cases, the other 
lengths taking "physical" values. That way the necessary number of points in the 
lattice will be reasonable (a few hundreds). Also, since the screening length is so 
large, the hard core of the macroions will not change the results drastically from the 
case of point particles so we will concentrate from now on on the Poisson-Boltzmann 
problem for point particles (Eq. f^q). 



III. RESULTS 



Equation (2.5) is solved numerically by an iterative method M. Using the Green 



function of the one-dimensional Laplacian 



G{z,z') = \\z-z'\ (3.1) 



the Poisson equation (2.1) can be written as 



(z) = -AttI I G{z,z')p{z')dz' (3.2) 







Starting with an arbitrary elec tric potential, one can compute the correspo ndin g 



density profiles using Eqs. ( p. 3D and derive a new electric potential using Eq. (^ 
then reiterate the process until a stationary solution is attained. In practice instead 
of using the new potential directly for the next iteration a mixing of the old and 
new densities is used. 



A. Generic results 

As stated before we have to consider very deionized systems in which the Bjerrum 
length I is smaller than 10 -6 A (the physical value of I for water at room temperature 
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is 7 A) . Figure [l] shows the density profiles of each species, the charge density, the 
electric potential and the electric field profile of a typical sample with the following 
parameters: I = 7-1CP 8 A, l g — 0.128 mm, h — 30 mm, Z — 100, a salt concentration 
Csait = 0.1 mMol/1 and a mean colloidal volume fraction ?7 co i = 0.12 (we consider 
that the particles have a hard core diameter a = 180 nm to express the colloidal 
density as a volume fraction in order to use units familiar with the experiments but 
we do not account for hard core effects in the Poisson-Boltzmann equation). 

The log plot of colloidal density profiles is similar to the experimental ones Q. 
In the bottom there is a slow decay whereas at high altitudes there is a faster 
barometric decay. Since we did not take into account the hard core of the particles 
in the theory we do not find the discontinuity in the density profiles near the bottom 
of the sample observed in the experiments due to the phase transition of the 
colloids from an amorph solid to a fluid. At very low altitudes (near the bottom, 
very high volume fractions) the Poisson-Boltzmann theory is not valid. 

The charge density profile confirms the results of Ref. || , that there is a strong 
accumulation of positive charges at the bottom of the container while there is a 
cloud of negative charge density floating at some altitude z* . There are clearly two 
neutral regions in the container: one at low altitude between the positive charge 
density at the bottom and the negative cloud, in which a non- vanishing electric field 
exists and a second neutral region at high altitude, over the negative cloud. The 
electric field in the lower region acts against gravity for the colloids therefore, as 
seen in the log plot of the colloids density profile, the decay is much slower than the 
one for an ideal neutral gas. Numeric results for other series of samples suggest that 
this electric field is proportional to Mg/Z. In the upper region the colloidal density 
drops exponentially as exp(— z/l g ) since the electric potential is almost constant 
and the electric field vanishes. 

Since the different densities vary with the altitude we can define a local screening 
length which depends on the altitude by 

X(z) = (4ttZ (Z 2 Pcol (z) + p+(z) + P -(z))y 1/2 (3.3) 

The two regions of the sediment are caracterized by a very different behavior of 
this local screening length. In the lower region the colloidal density is so high that 
Z 2 pcoi(z) ^> p+{z) + P-{z). In that region the colloids dominate the screening 
length. On the other hand in the upper region the colloidal density is very small 
and salt now controles the screening length which is then constant since at high 
altitudes the cations and anions densities are almost constant and equal to the 
mean salt concentration as seen in figure |l]. It is interesting to notice that electric 
charges accumulate in the intermediate region around z* where there is a change of 
regime, in agreement with macroscopic electrostatics principles. 

The preceding remark allows us to understand how the physical parameters (mean 
volume fraction of colloids, mass of the colloids, amount of added salt) will modify 
the altitude z* which separates the two regions. For example if we add more salt, z* 
will diminish since we reach sooner the regime where Z 2 p co \(z) < p + (z)+p-(z). We 
have computed the density profiles in several other cases changing the values of the 
parameters in order to find the depency of z* in these parameters. Our numerical 
results suggest that 



z*=-^^+a 2 ^3 (3.4) 

with ci = 0.15 ± 0.05 and a 2 = 1.0 ± 0.1. The preceding equation can be written 
in a more attracting way, introducing the screening length associated to the salt 
Asait = (47r/C sa it) -1 ^ 2 and the effective screening length associated to the colloids 
\f oy = (AnlZiN^/lg)- 1 / 2 , as 
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2* = A sa it y-ai + fl 2-^r^ (3-5) 

with ai — \[Akc\ = 0.5 ± 0.2. We do not consider here boundary effects: this 
equation is only valid if z* is smaller than h. The finite height h of the container 
will have the effect to "push" the negative cloud downwards if the parameters are 
such that z* approaches the top of the container. The same holds for the bottom 
of the container if z* is too small. 

Another quantity of interest is the size Az* of the negative cloud, defined as the 
mid-height width of the negative peak in the charge density profile (see figure |l|) . 
Since we know that at z* altitude, Z p co \{z*) is of the same order of magnitude as 
p+(z*) + p-(z*) — 2C sa it, the screening length at that altitude is proportional to 
A S ait- From basic electrostatics we know that the system will only tolerate charges 
over a length of order of the screening length, we deduce that Az* is proportional 
to A sa it- In fact the numerical results suggest also a linear dependency of Az* in l g : 

Az* = b 1 l g + b 2 \ salt (3-6) 

with bi = 5.0 ± 0.5, 62 = 0.7 ± 0.2, and the same restrictions concerning boundary 
effects as for the equation for z*. 



B. The apparent mass 



As we mentioned before, at high altitudes (larger than z*) the electric potential 
is almost constant and the electric field vanishes. From this it is clear that the 
colloidal density will decay as exp(— z/l g ) and there is no apparent reduced mass. 
Nevertheless let us notice that in the regime where the electric potential is almost 
constant in our calculations the corresponding colloidal volume fraction is smaller 
than 10 -9 . Such volume fractions cannot be measured experimentally. In practice 
the optical methods used in |Q allow to measure only volume fractions larger than 
10 -5 . A possible explanation to the apparent mass observed in the experiments is 
that for volume fractions higher than 10 -5 the asymptotic regime where the electric 
field vanish have not been reached yet: there is a residual electric field responsible 
of the observed reduced mass. 

In order to test this hypothesis we made a log plot of several colloidal volume 
fraction profiles restricting the plot to volume fractions higher than 10~ 5 (Figure |J). 
We computed the slope of the wing of the colloidal density to find an effective 
gravitational length I* which is higher than the actual gravitational length l g as 
observed in the experiments. Futhermore when we plot the colloidal volume fraction 
profile and the corresponding electric field profile together (Figure |J) we notice that 
for volume fractions higher than 10 -5 the electric field is not zero. 

The different plots in Figure [2] where obtained using different salt concentrations, 
so the sediment height (which is proportional to z*) varies. In this case we found 
that the apparent mass is a decreasing function of the height of the sediment, in 
agreement with the experiments. However the sediment height can be changed 
by changing other parameters like the mean colloidal density or their valency Z. 
Computing the apparent gravitational length I* as defined before for other series of 
samples, we found that the apparent gravitational length I* does not depend on Z 
or the mean colloidal density. In our model the ratio l* g /l g is only a function of the 
salt density. Figure (Q) shows the ratio l* g /l g as a function of salt screening length 

Asalt ■ 
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IV. COMPARISON WITH PREVIOUS APPROACHES 



As mentionned in the introduction the model presented above has motivated 
several studies both numerically [Q and analytically The purpose of this 

section is to compare our numerical results with the most achieved version of the 
theory presented in reference || . This theoretical approach is based on a constrained 
minimization of the free energy functional (|]^) assuming an exponential ansatz for 
the density profiles: 

jPcoi(z) = Nc ° l a exp(-az/l g ) (4.1a) 
'a 

P+ (z) = C sa it (4.1b) 

p-(z) = C sa it + ZNc ° l b exp(~bz/l g ) (4.1c) 

With this parametrization a — M* /M is the ratio of the reduced mass M* by the 
buoyant mass M of a colloidal particle, C sa i t denotes the fixed salt concentration, 
and N co \ is the fixed overall colloidal density per unit area, i.e. J^°° dzp co \(z) = 
N co \. The system considered in || is semi-infinite, z = corresponds to the bottom 
of the sample and h = +oo. a and b are the two variational dimcnsionlcss parameters 
of the theory, and the equilibrium density profiles p co \(z) and p~( z) c orrespond to 
the values of a and b that minimize the free energy functional ( p.2[ ). Following 
reference the minimization conditions express: 

- 1 ] (4.2a) 




v/l + (l-a)/ 7 



4 7 6 2 (a) 

-7+ p ±U k = ° 42b 

(a + b(a)Y 

where 7 = irZ 2 lN co \l g is the coupling parameter (I is the Bjerrum length introduced 
previously), and k = C sa \tl g /N co \ is the relative amount of added salt. Function / 
is defined by I(x) — J Q dy(ln(l + y))/y. 

Although equations ( |4.2| ) require a numerical treatment, it is possible to extract 
asymptotic expressions when the coupling parameter 7 is vanishingly small (strong 
gravitational coupling regime) or large compared to unity (strong Coulomb coupling 
regime) . Such an analysis is presented in reference ra] , and we only reproduce here 
the main features. When gravitational coupling is strong, 7 <C 1, the reduced 
mass is given by a ~ 1 — 37 (forall values of the salinity n) and therefore no 
mass reduction is observed in this regime (in agreement with the numerical results 
presented in reference [||). On the contrary, in strong Coulomb coupling regimes, 
7>1, quite a large mass reduction is predicted, even in low salinity regimes k <C 1 
(in such a situation the mass reduction is given by a ~ 1 + f ln 2 («(l + ^)) /{Z + 

1) + 0(1/7)). Our numerical results based on a free minimization of the functional 
( |2.2| ) show that it is indeed not the case, even though we observe nice exponential 
asymptotic behaviours at high altitudes. To emphasize this point we present in 
table I data obtained mostly in the strong Coulomb coupling regime. Excepted for 
the first value 7 = 7.2 10 -3 , which corresponds to the opposite strong gravitational 
coupling regime, the agreement between the theory presented in reference Q and 
the present free minimization is very poor. The effective mass predicted by the free 
minimization procedure corresponds to the actual mass within less than 0.1% in 
most situations, whereas the theory predicts effective masses that are only 2% of 
the actual mass! 
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To our opinion, the failure of the parametrization presented in reference |5j is not 
due to the exponential ansatz itself, but to the constraint of global charge neutrality 
applied to the asymptotic regime. The theory presented in reference [pi assumes 
that the profiles are exponential from the bottom of the sample to the top, as a 
result the free energy functional ( |2.2D has to be minimized with the constraint of 
global charge neutrality. However, the actual situation is quite different. If we 
refer to the experimental work done by Piazza et al. [Q, the exponential regime is 
reached only above a macroscopic layer of strongly interacting colloidal particles. 
Data presented in the previous section (see e.g. figure |l|) resulting from a free 
minimization of the functional also exhibit a dense macroscopic layer of colloidal 
particles in the bottom of the cell, and these profiles cannot be simply represented 
by a single exponential. This feature can be incorporated to the model suggested 
by Lowen, by splitting the cell into two parts. The upper part of the cell (above 
a given altitude "z " ) corresponds to the asymptotic region where the profiles can 
be accurately represented by an exponential, whereas below z the profiles are 
more complicated. As we can see, z Q is defined by the condition that the profiles 
are exponential above it. There is then no upper bound on the value of z Q and 
the asymptotic profiles should not depend on its precise value. As a result z 
can be chosen arbitrarily large. As a consequence, the part of the fluid located 
below z can be considered as a reservoir fixing the chemical potential of the ionic 
species /i co i and /i_ {p + is irrelevant since the local density p+(z) is held fixed, 
and is thus not a variational parameter). Although the full system must be charge 
neutral, the asymptotic part above z has no reason to be neutral. We are then 
lead to minimize the free energy of the upper part of the cell in the grand-canonical 
ensemble. Assuming that parametrization ( |4.lD is valid above z the minimization 
equation associated to the colloidal particles reads: 

dT\p co \,p + ,p-\ dN col 

^ = Mcol— k — (4.d) 

da oa 

where T[p co \, p+, p~] is now the free energy functional above z Q : 

/ + OC 
knTpiiz) [ln(A-p 2 (z)) - l] dz 

iEicol,+,-} -° 

"+oc ^ p-\-oc 



Mgzp co \{z)dz + - J ep(z)<5>(z) (4.4) 

and N co \ = J 2 + °° p co \(z)dz is the number of colloidal particles above z a per unit 
area. After some algebra, this minimization equation can be written on the form: 



(a-l){l + z: + zf}+a{f^-ln[X^'-^)}z, 



I Mcoi , / ,3 A'coia 



\k B T 



-z*/-, r, < -z*b/a a 1 z*(a 2 + b 2 ) 



= (4.5) 



where z* = az /l g . We c an easily check that when z* — we recover the first 



equation of condition (4.2). As z* can be chosen arbitrarily large in our model, we 



easily see that this equation implies a — 1 (no mass reduction) and: 

This new version of the theory is consistent with our numerical results predicting 
no mass reduction. 
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V. CONCLUSION 



A free minimization of the Poisson-Boltzmann theory used in references f2|,|5| 
have been performed in this article which lead us to conclude that this simple mean 
field theory does not predict any mass reduction contrarily to previous approximate 
minimization of the same functional. These new results arc fully consistent with the 
analytical results obtained in a two-dimensional case by Tellez ||. In particular, 
we observe the same condenser effect between the bottom of the container and the 
top of the dense region, resulting from a competition between electroneutrality and 
entropy of the microions. Data plotted in figure plgive a possible explanation for 
the experimental results obtained by Piazza et al. . Although in the asymptotic 
regime we observe no mass reduction, this regime is attained for very low values of 
the colloidal packing fractions, below the experimental resolution (10~ 5 ) in some 
situations. As a result, the residual electrostatic field can affect the profiles resulting 
in an apparent effective mass. 



9 



280 




Z (mm) z (mm) 

FIG. 1. From left to right, up to down, starting at the upper left corner, profiles of: 
counterions density, colloidal volume fraction in natural logarithm scale, coions density, 
charge density, dimensionless electric potential (j>, electric field E — —d<j)/dz. The param- 
eters used are: 1 = 7 ■ 1(T 8 A, l g = 0.128 mm, h = 30 mm, Z = 100, C sa it =0.1 mMol/1, 
fjcoi = 0.12 and a = 180 nm. 
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FIG. 2. Colloidal density profile in decimal logarithm scale for different salt concen- 
trations and restricted to volume fractions higher than 1CP 5 . Common parameters to all 



concentration in mMol/1 from left to right is 0.4, 0.3, 0.2, 0.1, 0.05, 0.04, 0.03, 0.02, 0.01, 
0.005. The apparent gravitational length l* g obtained from the slope of the low density 
wing is, from left to right in mm: 0.131 ±0.003, 0.131 ±0.003, 0.134 ±0.003, 0.140 ±0.004, 
0.151 ± 0.012, 0.155 ± 0.12, 0.162 ± 0.013, 0.172 ± 0.014, 0.187 ± 0.009, 0.217 ± 0.011. 
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FIG. 3. Colloidal volume fraction decimal logarithmic profile and the corresponding 
electric field profile in the case C sa it = 0.005 mMol/1, the other parameters being those of 
figure |i[ Notice that in the low density wing used to compute the apparent gravitational 
length lg the electric field is not zero. 
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FIG. 4. The ratio of apparent gravitational length by the actual gravitational length 
lg/lg versus the salt screening length A sa it, for two different values of the gravitational 
length. 
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Csalt 


Z 


M*/M : Theory § 


M* /M : free minimization 


7.2 


6 10 4 


3.27 10~ 4 Mol / 1 


100 


0.9915 


0.9994 ± 0.0001 


33.2 


0.33 


5 10~ b Mol / 1 


100 


0.0204 


0.986 ± 0.015 


33.2 


6.5 


10" 4 Mol / 1 


100 


0.0590 


0.999 ±0.001 


33.2 


325 


5 10~ 3 Mol / 1 


100 


0.3103 


1.00001 ±0.00001 


133 


6.5 


10~ 4 Mol / 1 


200 


0.0296 


0.9998 ± 0.0003 



TABLE I. Comparison between the reduced mass predicted by Lowen's theory and our 
numerical results from the free minimization. 
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